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ABSTRACT 

The study of the gas kinematic in the central 1.5 kpcxl.5 kpc region of M31 has revealed several surprises. The starting point of this 
investigation was the detection at the IRAM-30m telescope of molecular gas with very large line splittings up to 260kms~' within 
the beam (~40pc). In this region, which is known for its low gas content, we also detect an ionised gas outflow in the circumnuclear 
region (within 75 pc from the centre) extending to the whole area in X-ray. Relying on atomic, ionised, and molecular gas, we account 
for most observables with a scenario that assumes that a few hundreds Myr ago, M31 underwent a frontal collision with M32, 
which triggered some star-formation activity in the centre, and this collision explains the special configuration of M31 with two 
rings observed at 0.7 kpc and 10 kpc. The inner disc (whose rotation is detected in HI and ionised gas ([Nil])) has thus been tilted 
(inclination: 43deg, PA: 70deg) with respect to the main disc (inclination: 77deg, PA: 35deg). One of the CO velocity components is 
compatible with this inner disc, while the second one comes from a tilted ring-like material with 40deg inclination and PA=-35deg. 
The relic star formation estimated by previous works to have occurred more than 100 Myr ago could have been triggered by the 
collision and could be linked to the outflow detected in the ionised gas. Last, we demonstrate that the amplitude of the line splittings 
detected in CO centred on the systemic velocity with a relatively high spatial resolution (40 pc) cannot be accounted for by a possible 
weak bar that is roughly aligned along the minor axis. Although M31 has a triaxial bulge, there are no bar indicators in the gas 
component (photometry, no strong skewness of the isovelocities, etc.). 

Key words. (Galaxies:) Local Group, Galaxies: spiral, (Galaxies:) bulges, Methods: data analysis, Hydrodynamics, ISM: molecules 



1. Introduction 

The merging processes that probably took place during the as- 
sembly o f the galaxies of the Local Group is now actively stud- 
ied (e.g. iKlimentowski et all |201Q|) . While the Milky Way has 
not suffered any major merger for several billion ye ars and 
prese nts only a series of dwarf galaxy stellar streams (see lHelml 
2008), the other main spiral of the Local Group (M31) has had 
a more perturbed history. Giant stellar loops and tidal streams 
are observed in the surroundings, extending up to the neigh- 
bour _M322_wlrich_j^_mspected to have interacted with M31 
(e.g. Mc Connachie et all I2009I) . Many coherent structures ob- 
served around M31 are interpreted as the disruption of small 
dwarf galaxies, w hich are numerous in the vicinity of M31 (e.g. 
Ilbataetalll200l . 

At the same time, it has been known for a long time 
that the M 3 1 galaxy exhibits an unusual morphology (e.g. 
lArnl Il964t lHaas et all [l998l iHelferet all l200l . The young 
stellar population tracers i n the disc are all concentrated in 
the so-called 10 kpc ring dChemin et al 1 120091: Niete n et all 
2006; Azimlu et all 120111) . which is highly contr asted and su- 
perposed to only a few s piral structures (e.g. iNieten et ail 
2006} iGordon et all l2006h . One striking feature is the small 
amoun t of gas present in the cen tral region: neither iBraun et all 



(2009) nor IChemin et al .1 (120091) detect any significant HI com- 



ponent, while a sim ilar depletion is observed in CO-emission by 
INieten et all (120061) . It is usually described as a q uiescent galaxy 
with l ittle star formation SFR~0.4Mg vr" 1 (e.g. | Barm bv et al., 
2006; Tabatab aei & BerMiuiiserl^OldlAzTmlu et all 1201 ll) and 
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with an ultra-weak nuclear activity (del Bur go et all [2000). 
However, ion i sed gas is detected in the c entral field (e.g . 
[ Rubin & Fordl 119711: ICiardullo et all 119881; iBoulesteix et all 
ll987llBogdan & Gilfanovi 120081; iLiu et all 1201 Oh usually inter- 
preted in term of shocks, and lMelchior et al.T i 20001) has detected 
only a small amount of molecular gas ( 1 .5 x 1 4 M Q ) within 1.3' 
(305 pc in projection). There is n o obvious on-going star for- 
mation in this central reg ion (e.g. lOlsen et all l2006t iLi et all 
2009; lAzimlu et all 1201 1[) and the ionise d gas (approximately 
1500M o according to lJacobv et al.LlT985b can be accounted for 
by mass lost from evolving stars. 

In this "empty" centr al region, interest h as focused on the 
supermassive b lack hole dBacon et all 1 19941) . its activity (e.g. 
ILi et all 1201 lb and the circumnuclear region (inner few hun- 
dred parsecs) (ILi et all 120091) . iBarmbv et al.l (120061) observed 
the whole galaxy in the mid-infrared with the Infrared Array 
Camera (IRAC) on board the Spitzer S pace Telescope, rev ealing 
spectacular dust rings and spiral arms. iBlock et al.l (120061) stress 
the presence of an elongated inner ring with projected diameters 
1 .5 kpc by 1 kpc and propose a completely new interpretation for 
the morphology of this galaxy. Both rings at 1 kpc and 10 kpc 
are off-centred. Their respective radii do not correspond to what 
is expected from resonant rings in a barred spiral galaxy. The 
most likely scenario for the formati on of these rings is a head- 
on collision of the Cartwh eel typ e dStruck-Marcell & HigdonL 
[19931 iHorellou & Combes, 2001). Unlike the Cartwheel, where 
the companion is about l/3rd of the mass of the target (ma- 
jor merger), in Andromeda, the collision can be called a minor 
merg er, and produces m uch less contrasted rings in the main 
disc. IBlock et all (|2006) propose that the collision partner was 
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M32, with about l/10th of the mass (dark matter included) at 
the beginning. After stripping experienced in the collision, the 
M32 mass is now 1/23 that of the main target M31. The M32 
plunging head-on through the centre of M 3 1 has triggered the 
propagation of an annular wave, which is now identified with the 
lOkpc rin g, and a second wave propagates m ore slowly behind 
(see e.g. lAppleton & S truck-Marc elU. |l996), and would corre- 
spond to the inner ring. In addition, the inner ring has formed 
and is propagating in a tilted and warped disc, which accounts 
for its almost face-on appearance, in contrast with the inclined 
main disc of M 3 1 . This scenario explains why the cold gas has 
been expelled from the central region and also why there are 
shocks and hot gas. In this article we focus on the gas content of 
the inner ring and inside. 

We present CO(l-O) and CO(2-l) observations obtained with 
the IRAM-30m telescope in different positions located in the 
north-western part of the inner ring, and a few positions inside 
the ring. While the CO intensities are correlated in first instance 
with the Ab extinction, the velocity distributions are quite un- 
expected. In the inner ring the velocities are spread between 
-450 kms'and -150 kms _1 for a few positions while the sys- 
temic velocity (-310 kms~') is expected in this area along the 
minor axis. To better understand these velocities, we compared 
our measurements with the velocities available for the ionised 
and HI gas and found that they do not really match. While nei- 
ther the main HI warp of the disc nor a bar along the line of 
sight can explain the wide amplitude velocity splittings observed 
along the minor axis, the main configuration compatible with the 
data discussed here is to have a ring inclined with respect to the 
nuclear disc that would correspond in projection to the inner dust 
ring, which is detected in extinction and in infra-red emission. 
This ring could be attributed to a head-on collision with M32 
dBlock et all |2006) or to some accretion of gas from an M31- 
M33 gas bridge or tidal loop detected in HI dBraun & Thilkerl 
l200l . 

In Section [2] we describe the observations performed at the 
IRAM-30m telescope. In Section[3]we describe the reduction of 
the CO data. In Section [4] we present an analysis of the CO data 
and compare the molecular gas with the ionised gas, together 
with information from other wavelengths. In Section |5] we ex- 
plore various scenarios and propose one modelling that explains 
the observations. In Section [6] we discuss our results and con- 
clude. 

2. Observations 

The observations were carried out on 1999 June 13-15 and 2000 
July 14-17 with the IRAM 30-m telescope. Most of the obser- 
vations was made in the symmetrical wobbler switching mode, 
in which the secondary mirror nutates up to a maximum limit of 
+240 arcsec in azimuth. The beam throw was determined as a 
function of the hour angle in a way that t he OFF positions lay i n 
extinction-free regions, as described in iMelchior et al.l (|2000). 
It nevertheless occurred in two positions (M31D and M31G) 
that some signal was detected in the wobbler throw positions. 
As discussed in the Appendix [A] and indicated in Table [TJ this 
corresponds to some azimuth angles for M3 1G where one of the 
beam throws was close to the inner spiral arm in the north/north- 
western part of the field. For M31D the configuration could be 
more subtle and the OFF signal might have been caused by gas 
located on the far side, which was not detected in extinction. 
In the reduction process, the OFF signal is smeared out in the 
averaging process because equal weights are taken to reduce the 
observations to an equivalent beam of 24". It is still noticeable in 



M31G because the OFF signal was also present in the 2000 ob- 
servations. Near transit we had to use position-switching mode, 
taking an extinction-free OFF position located at a given posi- 
tion from the nucleus as indicated in the last column of Table 
[TJ We found the OFF signal for M31C for position-switch ob- 
servations, while no ON signal has been detected. These OFF 
detections are presented in Appendix lAl 

During the second epoch of observations small maps of 
seven points with a spacing of 5" were made for positions (A, 
C, D and G) to sample the CO(1-0) beam in CO(2-l) and to 
compute the CO(2-1)/CO(1-0) line ratio. Pointing and focus cal- 
ibration were regularly checked. In total, 10 ON positions in the 
inner disc of M 3 1 were observed in this way, as summarised in 
TableQJand presented in FigureQ] Table QJprovides (1) the name 
of the position observed, as used throughout the article, (2) the 
position angle in degrees of this position, (3) the offset in arcsec 
of the ON position with respect to the centre of M31, (4) the 
distance R of this position with respect to the centre of M 3 1 , (5) 
the right ascension, (6) the declination of this ON position, (7) 
the dates when the observations were made, (8) the total integra- 
tion time on the source T exp , (9) if a small map (seven pointings) 
was made for this position, (10) if some signal was detected in 
the OFF position, if yes, WSW (resp. PSW) indicates that some 
OFF signal was detected in the wobbler (resp. position) switch 
scans performed for this position, (11) the offsets in arcsec used 
for position-switch observations made near transit. 

We used four receivers simultaneously, two for I2 CO(1-0) at 
115GHz and two for 12 CO(2-l) at 230GHz. At 115 GHz each 
receiver was connected to two autocorrelator sub-bands (shifted 
by 40 MHz from each other) and each sub-band consisted of 225 
channels separated by 1 .25 MHz. At 230 GHz each receiver was 
connected to a filter-bank consisting of 512 channels of 1 MHz 
width. 



3. Data reduction 

The CLASSQ package was used for the data reduction. After 
checking the quality of each single spectra, the data were av- 
eraged with inverse variance weights. A first-order baseline 
was fitted to the resulting spectrum and subtracted. For a few 
spectra, a higher-order polynomial was carefully subtracted. 
Finally, the spectra were smoothed to a velocity resolution of 
3.2 km s 1 (resp. 2.6 kms 1 ) in CO(1-0) (resp. CO(2-l)). 

The maps are presented in Figure [2] and illustrate the ob- 
serving procedure. The spectra are located arbitrarily close to 
the beam position they correspond to. Spatial variations in the 
spectra agree relatively well with the extinction distribution. For 
M31G the CO intensity is clearly stronger in the area where the 
extinction is stronger. For M31A and M31D the intensity of the 
lines does not follow the intensity of the extinction exactly but 
we cannot exclude that some pointing errors and/or a clumpy 
distribution of the underlying gas explains the spatial variations 
of the CO distribution. 

The results of the fitting are presented in Table |2] A Gaussian 
function is fitted to each line to determine its area, central veloc- 
ity Vo, width <t, and peak temperature r pea k. The baseline RMS 
is provided for each line. We provide main-beam temperatures 
(unless specified otherwise) throughout this paper with B e ff = 
64.2+3 and F eff = 91+2 (resp. B eff = 42.±3 and F eff = 86+2) 
at 115 GHz (resp. 230 GHz). For the positions for which maps 
have been observed, the spectra were convolved to obtain an 24" 

1 Continuum and Line Analysis Single-dis Software, 
http://www.iram.fr/IRAMFR/GILDAS 
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Fig. 1. Right: Positions (J2000) observed in CO in the central part of M31. The circles i ndicate 1-mm (green ) and 3-mm (red) 
beams of the 10 positions observed superimposed on the As extinction map obtained in iMelchior et al.l (120001) . As discussed in 
this paper, the extinction with the central 2 arc sec is not well-de fined on this map. The (blue) cross indicates the centre of M3 1 
faf200o=00*42 m 44 J 371, 5nnm=41°16'08".34. ICrane et all (1 19921) /) Left: Velocity field measured in CO for the whole galaxy by 
Niete n et"a"D (|2006). A small rectangle in the central part indicates the field studied in this paper. This global view recalls the 
large-scale configuration of this galaxy: a steeply inclined disc at 77deg with a position angle of 35deg. 



beam ( Vl 1 2 + 21 2 ) for both lines, enabling a direct computation 
of the line ratio. 

The spectra are displayed in Figure [3] In the final spectra, 
the signal present in the OFF positions is mainly visible for the 
position M31G. It is at the systemic velocity and probably cor- 
responds to the gas in the main disc along the minor axis in 
the north of the M31G position (because the OFF positions are 
symmetric in azimuth with respect to the O N position). We use 
probably throughout this paper except for the lSaglia et al.l (l2010) 
data, whose velocities are formally in LSR, but the difference 
(~5kms _1 ) is negligible here. 

4. Analysis 

4. 1 . Characteristics of the CO emission and gas-dust 
connection 

The CO gas, which we searched for and discuss in this pa- 
per, is located mainly around the minor axis of the main disc 
of M3 1 . The detected molecular gas complexes correspond to 
the strongest Ab extinction complexes that we observed. For the 
positions without extinction (M31GB), we detected no signal, 
while no line was detected in the centre either. For M31E where 
the extinction is weak, we have a tentative detection at 4<x. The 
extinction wasc omputed assum i ng a fr action of foreground light 
of x = (see IMelchior et all [2000), 

which is probably cor- 
rect for the north-western part, where we observed, because an 



inclination of 45 deg i s usually assumed for the nuclear spiral 
dCiardullo et al.L[T988l) . We also computed the near-UV and far- 
UV extinction maps from archive GALEX data with the same 
method and found the same structures. In addition, we subtracted 
from the original GALEX images the bulge emission modelled 
in this way, and found no trace of UV emission in the inner ring. 



The position M31F, which exhibits a weak CO(1-0) signal 
with no detection in CO(2-l), corresponds to a position with 
low Ab extinction. The other positions are detected with high 
signal-to-noise ratios at least in CO(1-0). We computed the A^ 
extinction values for each position convolved with the differ- 
ent beams, as provided in Table [3] and there is no one-to-one 
correspondence with the intensity of the detected CO. Figure [1] 
also shows that M31D has a much stronger signal in CO than 
M31A while its A# extinction is weaker. One plausible expla- 
nation could be that the dust clumps corresponding to M31D do 
not lie in the same plane as M31A: if the foreground light (x) 
exceeds 0, we have underestimated the extinction. Assuming the 
Ico/Ab ratio measured for M31A applies to M31D, one would 
expect a peak intensity A r B = 1.56 for M31D, but it is mea- 
sured A™""- = 0.26. This configuration corresponds to a fraction 
of light in front of the dust x = 0.72, meaning that M31D lies on 
the back side of the bulge. It is also probable that the gas is very 
clumpy and that the non-linearity of the extinction somehow bi- 
ases our Ab estimate. 
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Fig. 2. Small maps performed at in CO(2-l) (left panels) and CO(1-0) (right panels) for the positions A, D and G. The green (resp. 
red) circles correspond to the FWHM of the CO(2-l) (resp. CO(1-0)) beams at the various positions observed. The corresponding 
offsets are: (0,0), (+5,0), (+2, +4) (-2, +4), (-5,0), (-2,-4), (+2, -4). The spectra are superimposed on the corresponding A^ maps. 
They are positioned arbitrarily close to the beam corresponding to the observations. 



Our detections are concentrated in an area0 of 415 pc x 
570 pc (in projection), located in the north-western part of the 
M 31 within 3.8 arcmin from the centre, corresponding to 880pc 

2 As in Melc hior et all ([2000), we assume a distance of M31 of 
780 kpc i.e. 1 arcsec=3.8 pc. 



(resp. 1.2 kpc if deprojected for a 45deg inclination). While the 
area explored is quite localised, the detected velocities span 
from -73 to -413 kms -1 , while velocities along the minor axis 
of a rotating disc are expected to be at the systemic velocity. 
M31K is the only position exhibiting this expected velocity, as 
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Fig. 3. Spectra with a detected signal. M31A, M31D and M31G (resp. M31E, M31F, M31I and M31J) spectra are reduced to a 243 
beam (resp. with the original beam) in CO(2-l) and CO(l-O). The y-axis provides the antenna temperature in K. 
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Table 1. Summary of observations 



Position 


PA 


Offset 


R 


RA 


DEC 


date 


X 


map 


Off 


Position 








(arcsec) 


(J2000) 


(J2000) 




(min) 




sig- 


switch 






















nal 


(arsec) 


M31A 


161.1 


-26,76 


80.3 


00:42:42.1 


41 


17:21.1 


13/06/99 
14- 

17/07/00 


630 
1160 


no 

yes 


no 
no 


122,321 
385,57; 
274,23; 
379,34 


M31C 


- 


0.0 





00:42:44.4 


41 


16:09.2 


2000 


256 


yes 


yes 


385,57 






















(PSW) 


M31D 


161.3 


-51,151 


159.4 


00:42:39.9 


41 


18:39.6 


14/06/99 


180 


no 


yes 


122,321 






















(WSW) 
















16/07/00 


112 


yes 


no 


404,-41 


M31E 


47.1 


72,67 


98.4 


00:42:50.8 


41 


17:15.0 


14/06/99 


360 


no 


no 


122,321 


M31F 


125.6 


-134,96 


164.8 


00:42:32.5 


41 


17:44.0 


15/06/99 


442 


no 


no 


122,321 


M31G 


121.8 


-171, 106 


201.2 


00:42:29.2 


41 


17:54.0 


15/06/99 


180 


no 


yes 


122,321 






















(WSW) 
















15/07/00 


175 


yes 


yes 


419,-7 






















(WSW) 


M31GB 


170.7 


-22,134 


135.8 


00:42:42.4 


41 


18:22.0 


15/06/99 


78 


no 


no 


122,321 


M31I 


123.8 


-172,115 


206.9 


00:42:29.1 


41 


18:3.6 
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41 
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Fig. 4. Superimposition of the G, I, and J spectra at the nominal spatial resolution, namely 21" for CO(1-0) and 11" for CO(2-l). 
The y-axis provides the antenna temperature in K. 



well as a component for M31G, M31I and M31J, and possi- 
bly M31D. Figure |4] shows a superimposition of the spectra ob- 
tained for M31G, M31I and M31J (and reduced to a 24" beam): 
the various velocity patterns appear in each spectra with dif- 
ferent relative intensities. Besides the velocity line detected at 
-145 kms -1 , a broadband signal is detected with velocities be- 
tween -450 and -250 kms -1 . As discussed below, this range of 
velocities cannot be explained by a single disc inclined at 45 deg 
in regular rotation, and they are the signature of a peculiar struc- 
ture. The line widths of the detected lines are quite different 
from one location to another, suggesting a spatial extension of 
the emitting areas. We assume a Galactic X C o = NhJIco = 



2.3 x 10 20 cnr 2 (K kms 1 ) 1 following lStrong et all (119881) . For 
the positions where we did not detect any signal, we provide a 
3 cr upper limit based on the dispersion (rms) computed on the 
baseline. Some CO(2-l) measurements (M31I, M31J, M31K) 
are at the limit of detection and are provided for a rough com- 
parison with CO(1-0) measurements. 

In the centre (< 2"), our Ab map does not measure the ex- 
tinction (owing to the method, which is based on ellipse fitting). 
The map of 8 /urn emission (after subtraction of the scaled stel- 
lar continuum at 3.6 fan) also displays a defect in the central 
part, so it is difficult to be c onclusive about the extinction in the 
inner arcseconds. However. iGarcia et al .1 (120001) . reiving on sim- 
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Table 2. Characteristics of the CO lines. CO(2-l) (resp. CO(2-l)) spectra are reduced to a 2.6 kms '(resp. 3.2 kms ') resolution. 



Position 


12 CO line 


I co (Kkms- 1 ) = 


V (km s -1 ) 


cr (km s 1 ) 


T eak 


baseline 


Nr 2 (cm -2 ) 






(beam size) 


f T,„ b dV 






(mK) 


rms 
(mK) 




(M G pc -2 


M31A 


1-0 (24") 


0.75±0.03 


-154.9+0.6 


28.7+1.3 


24.8 


2.1 


1.73xl0 20 


2.94 




2-1 (24") 


0.88±0.06 


-151.1+0.7 


22.9+1.9 


35.7 


4.7 


- 






1-0 (21") 


0.72±0.04 


-153.6+0.7 


28.7+1.7 


23.7 


2.4 


1.67xl0 20 


2.83 




2-1 (11") 


0.78±0.06 


-152.8+1.0 


24.9+2.6 


29.4 


5.5 


- 




M31C 


1-0 (24") 


- 


- 


- 


- 


3.7 


<2.55xl0 21 


<43.31 




2-1 (24") 


- 


- 


- 


— 


7.8 


- 




M31D 


1-0 (24") 


2.43±0.10 


-74.3+0.4 


23.7+1.2 


96.4 


7.2 


5.58xl0 20 


9.49 




2-1 (24") 


3.45±0.16 


-73.2+0.5 


19.7+1.1 


164.6 


15.5 


- 






1-0 (21") 


2.47±0.11 


-74.9+0.5 


23.6+1.4 


98.0 


8.7 


5.68xl0 20 


9.66 




2-1 (11") 


2.12±0.12 


-74.9+0.6 


21.1 + 1.4 


93.8 


10.6 


- 




M31E" 


1-0 (21") 


0.18±0.04 


-132.2+1.9 


15.8+3.1 


10.5 


4.1 


4.25xl0 19 


0.72 




2-1 (11") 


0.33±0.10 


-128.2+3.1 


15.3+6.9 


19.4 


13.1 


- 




M31F 


1-0 (21") 


1.56±0.14 


-360.6+3.3 


68.0+6.7 


20.6 


6.4 


3.59xl0 20 


6.11 




2-1 (11") 


- 


- 


— 


- 


16.3 


- 




M31G 


1-0 (24") 


2.98±0.18 


-147.4+0.4 


13.6+1.0 


205.9 


4.9 


6.86xl0 20 


11.66 




2-1 (24") 


2.39±0.14 


-146.7+0.5 


15.4+1.1 


145.2 


9.8 


— 






1-0 (21") 


3.01±0.23 


-146.7+0.6 


15.3+1.4 


185.2 


5.3 


6.92xl0 20 


11.77 




2-1 (11") 


2.18±0.20 


-146.6+0.9 


18.5+1.9 


110.6 


15.9 


— 




M31GB 


1-0 (21") 


- 


- 


- 


- 


7.7 


<5.29xl0 21 


<89.95 




2-1 (11") 


- 


- 


- 


- 


14.7 


- 




M31I 


1-0 (21") 


3.22±0.11 


-148.5+0.2 


9.3+0.4 


324.3 


12.5 


7.41xl0 20 


12.60 




1-0 (21") 


4.60±0.77 


-411.9+1.0 


34.1+3.1 


126.8 


12.5 


1.06xl0 21 


17.99 




1-0 (21") 


2.50+0.31 


-278.7+2.6 


46.3+5.4 


50.7 


12.5 


5.75xl0 20 


9.77 




1-0 (21") 


3.69±0.89 


-376.4+11.1 


88.5+16.5 


39.2 


12.5 


8.49xl0 20 


14.44 




2-1 (11") 


3.28±0.12 


-148.6+0.2 


7.7+0.4 


400.9 


29.0 


- 






2-1 (11") 


7.28±0.47 


-409.5+1.3 


40.0+3.5 


171.0 


29.0 


- 






2-1 (11") 


4.04±0.47 


-283.4+2.7 


43.5+5.0 


87.7 


29.0 


- 






2-1 (11") 


2.96±0.41 


-353.3+5.7 


46.6+4.5 


59.6 


29.0 






M31J 


1-0 (21") 






JJ. JZU.O 


361.4 


12.2 


2.95xl0 21 


50.08 




1-0 (21") 


1.06+0.11 


-145.5+0.6 


11.9+1.5 


83.8 


12.2 


2.45xl0 20 


4.16 




1-0 (21") 


7.34±0.57 


-350.1+7.5 


178.2+13.2 


38.8 


12.2 


1.69xl0 21 


28.70 




2-1 (11") 


15.22±0.78 


-412.7+0.9 


35.1+2.0 


407.0 


54.5 








2-1 (11") 


0.86±0.41 


-142.1+3.0 


9.0+5.2 


88.5 


54.5 








2-1 (11") 


1.94+0.69 


-268.6+7.1 


35.4+10.8 


51.4 


54.5 






M31K 


1-0 (21") 


2.66+0.14 


-300.8+1.1 


41.7+2.9 


59.8 


7.5 


6.11xl0 20 


10.38 




2-1 (11") 


2.61+0.43 


-292.5+4.9 


58.5+11.0 


42.0 


24.1 







Tentative detection. 



Table 3. A# extinction derived from the map based on the Ciard ullo et al.l ([l988) data as obtained in iMelchior et al.l d2000). For 
each configuration E(B-V) was then computed assuming Ry in the range [2.1,3]. These values were computed assuming a fraction 
of foreground light x — 0. The lco/E(B - V) ratio is then provided for the strongest lines. For M31J we also provide the ratio 
corresponding to the ring velocity. 



Position 


A B at 24" 


A B at 21" (11") 


E(B-V) at 24" 


E(B-V) at 21" 


E(B-V)at 11" 


WE(B-V) at 24" 


I C0 /E(B-V)at21" 




(mag beam -1 ) 


(mag beam -1 ) 


(mag beam -1 ) 


(mag beam -1 ) 


(mag beam -1 ) 


(K km s -1 mag -1 ) 


(K km s -1 mag -1 ) 


M31A 


0.095 


0.100(0.166) 


0.027+0.003 


0.029+0.004 


0.048+0.006 


28+5 


26+5 


M31D 


0.067 


0.076 (0.111) 


0.019+0.002 


0.022+0.003 


0.032+0.004 


129+22 


116+20 


M31F 


0.074 


0.081 (0.123) 


0.021+0.003 


0.023+0.003 


0.035+0.004 




69+15 


M31G 


0.182 


0.193 (0.256) 


0.052+0.007 


0.055+0.007 


0.073+0.009 


59+11 


56+11 


M31I 


0.210 


0.219 (0.245) 


0.060+0.008 


0.063+0.008 


0.070+0.009 




52+ 8 


M31J 


0.203 


0.208 (0.241) 


0.058+0.007 


0.060+0.008 


0.069+0.009 




219+33(18+4) 


M31K 


0.138 


0.141 (0.159) 


0.040+0.005 


0.040+0.005 


0.046+0.006 




67+12 



pie modelling of Chandra data, have estimated Ay = 1.5 + 0.6 
in the central 3-arcsec region. Our millimetre observations ex- 
clude the presence of CO(1-0) (resp. CO(2-l)) at the 1 cr level 
of 3.7 mK (resp. 7.8 mK). On the basis of this upper limit, the 
gas present within 80 pc from M31's centre does not exceed 



43M /pc -2 (3 cr), assuming a conservative value for Xco typi- 
cal of the Galactic disc. This is an extremely small amount com- 
pared toJhej*^masscomputed in the central region of the Milky 
Wav (IOkaet allll998h . 
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Table 4. Line ratios of the complexes reduced to a 24" beam. We 
provide both the ratios of the line intensities Tu and the ratios of 
the peak temperatures R. 



Complexes 


rn 




R 


M31A 
M31D 
M31G 


1.17= 
1.43; 
0.80d 


.0.40 
c 0.19 
11.45 

~n 

"0.13 


1.44±" 
1.71±°- 3 5 

07i+8:li 

U -' ±:E 0.05 



We used the three positions with observations reduced to a 
24" beam to compute the CO(2-l) to CO(l-O) line ratio. As dis- 
played in Table |4] we computed both the ratio of the integrated 
intensities and the ratio of the peak temperatures, which are com- 
patible within error bars. M31G exhibits a line ratio compara- 
ble with the mean valu e found for the Milky Way spiral arms 
dSakamoto et all 1 1997k 

M31A and M31D exhibits CO line ratios higher than 1. 
Given the N# 2 column densities computed (respectively 5.6 x 
10 20 cirT 2 and 1.7 x 10 20 cm -2 ), the gas is probably not optically 
thin. In addition, the velocity dispersion of CO(2-l) is systemat- 
ically smaller than that of CO(l-O), suggesting that these com- 
plexes are probably externally heated and composed of dense 
gas. This hypothesis is compatible, as discussed below, with the 
combined pres ence of an intense st ellar radiation field caused 
by bulge stars (Stephens et al., 2003), the outflow detected in X- 
ray in this area dBogdan & Gilfano v, 2008), and the ionised gas 
dJacobv et all 11985 ). This supports the view that no star forma- 
tion is currently happening inside the inner ring of M3 1 . 

4.2. Ionised gas 

Th e ionised gas in t he ce ntral regions of M 3 1 has been mapped 
by ICiardullo et al.l d!988l) . who produced the highest S/N ratio 
map of Ha+[NII] based on the compilation of a five-year nova 
survey published so far. It exhibits, as displayed in the top left 
panel of Figure|5] a turbulent spiral, which appears more face-on 
than does the main disc. These authors described it as lying in a 
disc, warped in a way that the gas south of the nucleus is viewed 
closer to face-o n than the gas in th e northern half of the galaxy. 
As discussed in lBlock etatl {2006), this arc feature corresponds 
to the sout h-west part of the of fset ring detected with Spitzer- 
IRAC data. lRubin & Fordid 197 lb measured that in this area [Nil] 
is three times stronger than Ha, which supports the shock hy- 
pothesis as also discussed by Jacobv et al. ( 1985). According to 
iLiu et alj d2010l) . half of the kinematic temperature of hot gas 
in the central bulge is accounted for by the stellar dispersion, 
but additional heating is expected from type-la supernovae, even 
though the question of the iron enrichment i s not entirely clear . 
A head-on collision with M32 suggested bv lBlock et alJ 12006) 
could account for the additional heating and at least contribute 
to it. 

In parallel, Boulesteix et al. (1987) observed this same area 
with a Fabry-Perot spectrograph to map the velocity field, and 
produced an [Nil] map, which is xdisplayed in the top right 
panel of Figure|5] The centre has been masked for technical rea- 
so ns. Besides the lower r esolution and lower S/N ratio than for 
the lCiardullo et al.l d l988) map (because of a smaller integration 
time), it is instructive to observe that it exhibits different pat- 
terns compared with the Ha+[NII] map. It is striking that the fil- 
ament, which crosses the position M31A, is not clearly detected 
in [Nil], so it should be mainly associated to Ha component. 
The comparison of these two maps enhances the fact that ob- 



viously forbidden lines and Balmer lines do not sample exactly 
the same regions: (1) the gas does not have the same Ha/[NII] 
ratio everywhere; (2) the extinction additionally complicates the 
analysis. 

4.2.1 . Velocity field of the central bulge (1 .5 kpcxl .5 kpc) 

To compare the ionised gas with our molecular detections, we 
tried to reconstruct the ionised gas velocity field from data pub- 
lished in the literature. Various lines have been used, but only 
a single velocity has been measured for each spectrum. In the 
bottom right panel of F igure |5] we display the velocity map of 
Boulestei x et al.l d!987l) . The best S/N ratio is of course observed 
where the [Nil] intensity is strongest (cf. top right panel of Fig. 
0. It exhibits an irregular disc in rotation, with obvious pertur- 
bations along the minor a xis. Planetary nebulae are numerous in 
this field (537 detected bv lMerrett et all I20061) and add noise to 
this velocity field. As discussed below (see Tabled, we estimate 
here a position angle of 70deg (-20de g for the kinemati c mino r 
axis) compared to 40deg claimed by iBoulesteix et alJ dl987l) . 
ISaglia et all ([2010) interpreted these perturbations as counter- 
rotations. In the bottom left panel of Figure [5] we display the 
slit measurements we gathered from the liter ature, using Dextep 
when necessary, nam ely H/? and TOm l fromlSaglia et all yOlO); 
[On] and [Neml fromlCiardullo et al l (119881) : Ha and [Nn] from 
lRubin&For"dl(fl97ll) . We also superimpose our CO velocities, 
and find a very marginal agreement. In Figure |6] we superim- 
pose the two figures. We find a good overall agreement for the 
various ionised gas measurements with complicated features in 
the central part and along the minor axis, with the most striking 
discrepancy at position angles of 128 deg and 142 deg. 

It is difficult to understand the discrepancies observed in the 
ionised gas velocity field, as the techniques are quite different 
and the IBoulesteix et all dl987l) data suffers relatively low reso- 
lution (even though excellent given the fact that it has been ob- 
tained in 1985!). We can consider several explanations: (1) the 
various lines might have different relative ratios from one region 
to another (as seen in the top panels of Figure [5} and might be 
affected differently by extinction (see Figure [Q. (2) It is possi- 
ble that the ionised gas, like the molecular gas, exhibits several 
velocity components, a nd it has been not explored so far, but by 
IBoulesteix et al.l ([1987), who mentioned line splittings greater 
than 30 kms 'in the central region. 

In addition, the CO velocities do not follow the regular pat- 
tern. M3 1 A and M3 ID do not really match the ionised gas veloc- 
ities. However, some CO positions might have a velocity (or at 
least one component) compatible with the ionised gas, namely 
M31E, M31F, M31K, M31I, M31J and M31G. This suggests 
that only part of the molecular gas is kinetically decoupled from 
the ionised gas. 

Also, we wanted to figure out whether similar li ne splittings 
were present in the ionised data. As discussed by Sagli a et all 
d2010h . the intrinsic velocity dispersions of the gas is smaller 
than 80 kms _I , while the in strumental resolution ac hieved by 
these authors is 57 kms" 1 . [Boulestei x et alJ dl987l) . who de- 
tected line splittings larger than 30 kms~' in the central area, 
reached a spectral resolution of 14 km s _I . It is thus challenging 
to detect line splittings of this order of magnitude, but such split- 
tings might be underlying and explain part of the discrepancies. 
As we observe large (> 200 kms -1 ) line splittings in CO, we try 
to investigate if such double components can be detected from 
the ionised gas. 

3 http ://dc . zah.uni-heidelberg.de/sdexter 
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Fig. 5. Ha+[NII] (left) and [Nil] (right) intensities from ICiardullo et all d!988l) and iBoulesteix etalj d!987l) are displayed with an 
arbitrary normalisation in the top panels. The bottom left panel present s the CO velocities measured in this paper, superimposed on 
the gas kinema tics based on (1) H/3 and [Om] lines measured in slits by Saglia et al. j|2010l) . (2) [On] and [Nem] measured in slits by 
ICiardullo et al.l d 1 988b and (3) Ha and [Nn] measured in slits by Rubin & Fordldl971l) . The bottom rig ht panel displays t he ve locity 
field obtained with [Nn] /16584A observations (cf top right panel) based on a Fabry-Perot device by IBoulesteix et ali (|1987). On 
each panel, the positions of our CO observations are indicated with circles corresponding to the beams. 



Table 5. Line splittings detected on [Om] and H/3 lines from Saglia et al.l (I2010I) . We provide the velocities, signal-to-noise ratio of 
the integrated line and [O m]/H6 line ra t io (log u 
single values published by Saglia et al. (2010). 



the integrated line and \O m~\fH0 line ra t io (lo gm) for each detected component (1 and 2). We also give for comparison purposes the 

~| (Hi 



offsets 


R(") 


PA(deg) 


v [OUi]5007 (kmirl) 


S/Nj 


[Qm]/H0 




S/N 2 


[Om]/H£ 


v « as 

Saglia 


-15.2,16.9 


-22.8 


138 


-548.6±12.0 


4.7 


51+0.20 


-360.6+16.7 


5.2 


-0 22 +u lu 

-0,12 

0.79+f 
0.863$* 


-469.1+8.4 


-13.5,-12.2 


-18.2 


48 


-462.8+6.0 


13.1 


-Q.13 
n 9Q+0.09 

-Q-Q8 


-234.0+4.0 


15.4 


-323.7+7.4 


-12.0,-10.8 


-16.2 


48 


-449.5±9.1 


10.4 


73 + ° - 8 

—0.16 


-225.4+5.2 


12.2 


-314.8+5.4 


-9.3,-8.4 


-12.5 


48 


-485.8+9.9 


4.9 


-375.0+6.7 


3.9 


-0 23 +0 - 20 


-418.5+10.7 


-8.1,-7.3 


-10.9 


48 


-454.3+14.4 


3.3 


15 +023 
-0.25 


-272.5+13.3 


5.3 


041 +0.22 


-314.4+9.2 


-1.9,-5.9 


-6.2 


18 


-319.9+20.0 


3.3 


0.141 | 


-144.8+15.9 


3.6 


19+ - 32 

"■"0.28 


-342.7 


-0.2,0.9 


-0.9 


168 


-281.3+8.9 


8.7 


0.78!]| 


-85.9+16.2 


3.2 


o oo +0 - 72 

"'""0.41 


-241.8+10.0 


0.1,-0.7 


0.7 


168 


-319.9+19.8 


3.3 


-0 20 +o f4 

-0J8 

60 +0 22 
-0 06 +o ft 


-144.8+16.1 


3.7 


-0 04 + ° 20 

1 04 +lnT 
I -"^-0.40 


-221.7+11.6 


0.7,0.6 


1.0 


48 


-332.6+11.0 


6.8 


-64.9+7.5 


7.5 


-313.6+7.8 


1.1,1.0 


1.5 


48 


-322.3+13.0 


4.9 


-56.0+3.5 


18.2 


92+ 032 


-330.6+7.8 


1.2,0.3 


1.2 


78 


-321.5+13.2 


5.2 


20 +0 26 


-78.7+4.8 


12.7 


0-57 +° 7 22 


-327.8+11.9 


8.8,-2.9 


9.3 


108 


-598.0+4.7 


14.1 


92 +0 37 


-319.2+14.0 


8.1 


0.18!»! 2 


-343.6+13.9 


11.2,-3.6 


11.7 


108 


-593.7+5.5 


8.4 


i 17+0.56 
l.l / „,„ 

-0 16 +0 ' 7 


-373.3+29.2 


3.7 


0.28!»;|I 


-376.8+18.1 


5.3,16.4 


17.2 


18 


-575.1+3.8 


6.4 


-444.8+35.2 


4.5 


-291.4+6.6 
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Fig. 6. 6.7 x6. 7 ionised gas velocity field centred on M31. 
Superimpo sition of the various gas velocities (see Fig.[5]bottom 
left) on the Boulest eix et alJ (Il987t) Fabry-Perot-derived velocity 
map. 



41°16 , 10" - 



41*1 6W - 




[ig h 4Z m 46'0 45'0 44?D 43?Q 

Fig. 7. 40 x40 ionised gas velocity field centred on M3 1 . Line 
splittings detected in [Om] 50007 A are superimposed on th e ve- 
locity field measured on ionised gas by Sa glia et all (12010 ) and 
bv ldel Burgo et al.l(l2000l) . This velocity field is quite unexpected 
and do not display any clear rotation pattern. 



Optical spectra of M31 are heavily dominated by the stel- 
lar continuum of its bulge. Therefore, in order to study emis- 
sion lines in that spectral region we approximated the calibrated 
sky subtracted spectra of M 31 kindly provided by R. Saglia by 
stellar pop ulation models using the N Bursts full spectral fitting 
technique (Chiling arian et al.ll2007b ;al. We excluded narrow re- 
gions around H/3, [Om], and [Ni] lines from the fitting. Then 
we subtracted the best-fitting stellar population models from the 
original spectra and analysed the fitting residuals. This allowed 
us to reliably measure parameters of these rather faint emission 
lines which are barely visible in the original data. We then in- 
spect visually the spectra to determine the positions where high 
S/N line splittings is present. We then fit a two Gaussian func- 
tions with the CLASS package on the [Om]5007A line. Template 
mismatch and the relative low S/N affect the measurement of the 
H/3, which is weaker and not always detected. For none of our 
measurements (but one), the [NI] line is detected. Our detection 



are summarised in Table [5] and Figure [7] We identify 14 posi- 
tions with a double component detected in [Om] and also in H/3, 
within 23" (~86pc) from the centre, but we do not detect any 
line splittings close to our CO observations or at a similar radial 
distance. 

4.2.2. Velocity field of the circumnuclear region 
(150pcx150pc) 

Figure [7] displays the velocity field measured in ionised gas 
of the circumnuc lear region. This includes slit data from 
Sagli a et al.l (|2010|) as well as two-dimensional spectroscopy of 
Idel Burgo et al.l d2000l) for the measurements of single veloc- 
ity components, a nd our line splitting measurements based on 
Saglia etal. (2010) data. It is striking that unlike Figure[6]this ve- 
locity field does not exhibit any clear rotation pattern in the cen- 
tral field. While the centre presents a spot at the systemic veloc- 
ity, the whole area has velocities smaller than the systemic veloc- 
ity and is approaching the observer. This coherent flow of ionised 
gas is decoupled from the stellar kinematics (e. g . iBender et all 
120051; Idel Burgo etafl 120001; ISaglia et all 12010)) , and resembles 
an ionised gas outflow. This is confirmed by the fact that the 
two indepen dent data sets are in good agreement. We should 
stress that the Saglia et all (1201 Oh data are characterised by an in- 
strumental resolution of 57 km s _1 , and their systemic velocity is 
shifted by +23 km s _1 with respect to the other studies. However, 
we checked that this cannot account for the systematic blue shift 
observed with respect to the centre. The mechanisms that rule 
the feed-back from the AGN are still largely unknown and in the 
case of M31, we are only tracing a relic activity of its known 
black hole. Considering that M 3 l's black hole is not active, two 
main mechanisms can be considered to explain this outflow. (1) 
It could be directly associated to the past stellar activity detected 
in the central regions: a 200 Myr old ste l lar dis c with a 24OOM 
mass has been detected bv lBender et alJ (120051) within a fra ction 
of arcsec from the centre. Inside a 2"disc. lSaglia et all (1201 Ol) has 
also detected a stellar component younger than 600 Myr corre- 
sponding to a mass smaller than 2xlO 6 M : it could be compati- 
ble with a 100 Myr component with a mass of 1O 6 M . However, 
we lack of constraints to conclude for a star formation origin 
with a typical dynamic time r^,.,, = size / velocity of this out- 
flow. The typical velocity of an outflow (e.g. lRupke et all [20051) 
is about 400 kms~', a Td yn = 100 Myr old starburst would then 
correspond to a size of the emitting area of 40kpc. It might cor- 
respon d to the X-ray outflow detected by Bogdan & Gilfanovl 
(2008) (see Figure [TT] lef t panel). (2) As d iscussed by various 
authors (e.g. [David et all l2006l Hoi 120091) . the ionised gas in 
galaxy bulges can be accounted for by mass loss of evolved 
stars and SNIa could play a key role. The prompt SNIa compo- 
ne nt (which represents 50 % of the SNIa population) discussed 
by iMannucci et al.l d2006l) could hence be indirectly linked to 
the 100-200Myr star formation episodes quoted above and con- 
tribute to this outflow. 

As discussed in the following, a possible collision with M32 
can explain the ring structures of M3 1 where the CO has been 
detected, and such an event could have triggered a star formation 
activity in the central part and support the SNIa explanation for 
the outflow. In addition, this inner outflow could be connected 
with the outflow detected in X-ray by (Bogdan & Gilfanovl 
2008) on a larger scale and discussed in Sect. 14.41 

Several (8/14) of the positions with line splittings have a 
component with a large [Oin]/Hj8 ratio. Such a large ionisation 
is compatible with planetary nebulae, del Burgo et al. (2000) 
observed several high-ionisation "clouds" in this area. These 
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Fig. 8. Superposition of the various gas velocities (see Fig. [5] bottom left) (left:) on a simple model of a galactic disc with an 
inclination of 45deg, a position angle of 40deg and a systemic velocity of -310 kms -1 ; (right:) similar disc with variable position 
angles 




Fig. 9. Image of the HI emission in the very centre, avoiding 
+ 100km/s around the systemic velocity of -310km/s. This cut 
has been made to subtract the large-scale HI emission seen in 
projection superposed to the centre, expected near the systemic 
velocity. The contours are the dust emission from Spitzer-IRAC 
(8/im image where a scaled version of the 3.6 fim image has 
been subtracted to remove the stellar photospheres emission) 
( IB lock et all |2006|). The HI cube is from Braun et al (2009). 



authors discussed that the intensity of 3 of these sources is 
much larger than thos e of planetary nebulae detected in M3 1 by 
ICiardullo et al.l (120021) . However, in this central region, it might 
be possible that these sour ces are multiple. Intere stingly, one of 
the sources (D) detected bv ldel Burgo et al.l (12000) exhibits a line 
splitting (-270,-527)kms~' with an amplitude comparable with 
ours but very close to the centre, even though none of the veloc- 
ities really match. Also, the strongest component we detect in 
the inner ar csec region matches in first approximation with the 
source A of Idel Burgo et al.l (120001) . 

A comparison with SAURON data (M. Sarzi, private com- 
munication) taken in this region of M3 1 has shown that our mul- 
tiple lines detections are compatible with planetary nebulae. The 
slit most probably biases the line ratio determination, as different 
parts of the PSF are sampled. 




Fig. 10. Velocity field of the HI emission, corresponding to 
Figure [9] avoiding large-scale HI emission. 



4.2.3. Modelling of a tilted inclined disc 

We create a velocity field map with the CCD VEL task 
within the NEMO so ftware package (Teuben U 19951) . Following 
Ciardullo et al. (1988), we consider a tilted ri ng model with an 
inclin ation of 45deg, a position angle of 40deg dBoulesteix et all 
1987) and a systemic velocity of -3 10 km s _1 . In order to under- 
stand the behaviour of the gas in the central part of M31, we 
superimpose the velocities measured in this paper (ionised gas 
and CO gas) on this modelling in Figure [8] (left panel). Even 
though there is a veloc ity field, as shown with ionised gas by 
iBoulesteix et al.l dl987l) and displayed in the bottom right panel 
of Figure [6] it is not regular and does not exhibit a well-defined 
zero- velocity curve. In addition, the velocities of the CO mea- 
surements do not fit with such a regular speed pattern. Varying 
the velocity angle in the centre, as displayed in the right panel 
of Figure [8] mimics the velocity pattern in the inner region 
(< 10 arcsec) but does not explain the minor axis configuration. 
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4.3. HI emission 

The detailed maps from lBraun et al.l d200 9) have revealed clearly 
the HI deficiency in the central parts of M3 1 . There is, however, 
some HI emission in the central kiloparsec, but most of it comes 
near the systemic velocity, and is likely to be projected emission 
from the external tilted gas o rbits, which are warp ed to an al- 
most edge-on inclination (e.g. Corbelr Tet al.L 120 lOl) . To subtract 
this projected large-radii emission, and better see the residual 
coming from the actual centre, with a large velocity gradient, 
we have summed all channels avoiding ±100km/s around the 
systemic velocity of -310km/s. A weak signal can then be dis- 
tinguished, with a morphology of an incomplete ring (Figure |5), 
with emission corresponding to the east side of the inner ring 
delineated by the dust emission revealed in the 8/mi Spitzer- 
IRAC map. Figure [9] shows the dust contours superposed on 
the residual high-velocity central HI emission. In this picture, 
we can see that the main HI residual component still follows 
the large-scale (NE-SW) arms seen in projection, and coincid- 
ing with the dust features, in the direction of the major axis of 
M31. However, there remains a weak component perpendicu- 
lar (NW-SE) to it. The main concentration of this residual HI 
emission is well aligned with the north-east part of the dus t ring 
and it could correspond to the weak component seen by Brinks 
(Il983h on the minor axis position velocity diagram (his figure 
lb). However, the remaining parts of the ring is devoid of HI 
emission; the atomic gas must have been transformed into the 
molecular phase in the dense parts of the ring. As discussed 
in Appendix iBl Figure IbTI reveals that most of the CO strong 
emission at high velocity in the centre has no HI counterpart. 

The HI component associated to the dust ring is, however, 
too weak to be seen in the velocity map, as shown from Figure 
[TUl The velocity field in the HI selected with high speed with re- 
spect to the systemic velocity, is still compatible to the "normal" 
main disc, and similar to the velocity field found with the ionised 
gas in Ha, although with possible perturbations in the vicinity of 
the ring. It is difficult to determine an accurate position angle of 
the HI isovelocity map (Figure IT0T > to discriminate between the 
nuclear disc and the main disc, but the intensity map (Figure |9]l 
does not really favour a nuclear disc component. 

4.4. Other wavelengths information 

While the centre of M3 1 hosts a massive black hole with a mas s 
of 0.7 - 1.4 x 10 8 M o dBacon et all 120011: iBender et all 120051) . 
it is one of the mo st underluminous supermassive black hole 
(Garci a et al ., 2010). The A-star clust er, detected in the third 
component (P3) of M 31's nucleus by Ben der et alj ([2005), can 
be associated to a recent star formation episode, typically a 
single burst, which occurred 200 Myr ago, with a total mass 
in the range 10 4 - 1O 6 M , corresponding to an accretion rate 
of 10~ 4 - 10~ 2 M o yr _1 . It might be triggered by the possible 
frontal collision with M32 (Blo ck et aiTT2006). The detection 
of an ioni sed gas outflow in X-rays a long the minor axis of the 
galaxy by Bog dan & Gilfanovl (120081) . perpendicular to the main 
disc, could be linked to this recent star formation activity and to 
the possible outflow detected in the optical (Figure|7]i. However, 
this is still controversial as such a burst d oes not have enough 
energ y to power the galactic wind required dBogdan & Gilfanovl 
2008). In the left panel of Figure[Ti] Ha+[Nn] contours are su- 
perimposed (with the same resolution) on the C handra soft X- 
ray emission map from lBogdan & Gilfanovl (120081) . The relative 
intensity of the outflow on both sides is compatible with the in- 
tensity of Ab extinction: the NW side is more extinguished than 



the SE side. This is compatible with the modelling described in 
the next section. Last, the right panel of Figure Qj] displays su- 
perimposition of the contours of the A# extinction map on the 
PAH-dust emission at 8 //m detected by Spitzer. The extinction 
features match exactly the ring detected at 8 yum. 

5. Interpretation 

The molecular, atomic and ionised gas exhibit diffe rent radial 
distrib utions in disc galaxies as first discussed by iKennicutH 
(119891) (see also lBigiel et~aTll2008l for a recent review). The HI 
gas is known to extend at muc h larger radius. M3 1 has HI gas ex- 
tendin g at least up to 40 kpc dCorbelli et alll2010tlChemin et al. , 
2009) with high velocity clouds up to 50 kpc dWestmeier et al. , 
2008). As discussed below, the projected warp dominates the HI 
emission in the central part (see also|B]i. The CO and ionised gas 
are usually more conc entrated. In M3 1 , w hile CO is known to be 
depleted in its centre (Nieten et al., 2006), the ionised gas is de- 
tected in the inner parts (Ciardullo et al.lll988HBoulesteix et all 
119871) and is strongest there dDevereux et al.l 11994 . 

While the atomic and ionised gas exhibit the presence a per- 
turbed disc, the molecular gas detected in CO displays unex- 
pected kinematic signatures, with significant line splits close to 
the minor axis and a very weak when present signal close to the 
systemic velocity. Several scenarios could be invoked to explain 
the existence of two well separated high S/N velocity compo- 
nents ( A V = 260 km s~ 1 on a scale of 40 pc) traced by the molec- 
ular component (CO emission) in the centre of M3 1 . One com- 
ponent is in the sense of the expected rotation, the other com- 
ponent is counter-rotating. The widths of these two components 
are also very different, as displayed in Figures|3]and|4] 

The ionised gas also exhibits unexpected kinematic features 
with a face-on outflow in the circumnuclear region (detected in 
the optical), extending (in Xray) to the 1.5 kpc x 1.5 kpc region 
studied in this paper. 

In the following, we explore several possible scenarios to 
explain the observations. In Sect. 15.11 we discuss the expected 
effects of a large scale warp. In Sect. 15.21 we summarise the ex- 
pected signatures of a bar and demonstrate that they fail to match 
the CO observations. In Sect. 15.31 we remi nd the scenario of a 
head-on collision with M 32 as proposed by iBlock et al.l (12006) 
and show with a simple modelling that it can account for the 
various observables. 

5. 1 . Scenario 1 : large scale structures warp and tidal 
streams 

We discuss the arguments that the second velocity component 
detected in CO could be due to the superimposition of gas at 
large scale onto the central regions. These possible configura- 
tions are not related to the outflows detected there in the ionised 
gas. 

5.1 .1 . Large scale warp 

We could think that the second component, the counter-rotating 
one, is just observed in projection, and comes from the exter- 
nal warp observed in the outer parts of the M31 disc, with a 
different inclination and position angle than the normal M31 
disc. The warp is vis ually conspicuous, and remarkable in HI 
dCorbelli etall I2 010). Far t her in the north-east disc (beyond 
10 kpc). ICasoli & Combesl (Q~988) already noted that CO emis- 
sion could come from the main disc and the warped component, 



12 



A.-L. Melchior and F. Combes: Molecular gas in the inner 0.7 kpc-radius ring of M31 




Fig. 11. Left: S uperposition of the Hq+fN nl map (upper left panel of figure [5} of ICiardullo et al.l (1 19881) on the Chandra soft X-ray 
emission from Bog dan & Gilfanovl (120081 see their Figure 7). The Ha+[Nii] map has been smoothed with the same smoothing 
lengths as those used for the Chandra map (A. Bogdan, private communication). Righ t: Superpositi on of the As extinction on the 
8/im Spitzer-IRAC map after subtraction of a scaled version of the 3.6 fim image ( iBlock etail 12006). The Ab contours are fixed to 
0.05 (blue) and 0.15 (white). 



both on the same line of sight. In that case, the CO emission 
(T* ~ 40mK) was coming from material distant by more than 
16 kpc from the centre. 

This scenario is, however, impossible, because (1) the ex- 
pected molecular gas has negligible emission at the large radii 
where this warped component should be (e.g. iNeininger et all 
Il998l) : (2) the corresponding projected velocity on the line of 
sight close to the centre is expected around the systemic velocity. 
It is likely that the gas at long distance (16-30 kpc) is on nearly 
circular orbits and has no radial velocity when projected to the 
centre. Accordingly, the corresponding warped material was al- 
ways found at systemic velocity. Whatever inclination it has, the 
gradient in velocity across a small region of 40pc should be neg- 
ligible. Ad hoc hypotheses with HI infall/outflow with important 
radial velocity are then required, which are incompatible with 
the HI observations (e.g. Corbelli et al 2010). We therefore think 
that, unless the gas is completely out of e quilibrium (which is 
not r eally sup ported by t he observations of Corbelli et all 2010; 
ICasoli & Combesl Q988), the peculiar component cannot come 
from the outer warp. However, we have some tentative detections 
close to the systemic velocity (M31K, M31I, M31G, M31J), 
which could be accounted for by the external warp. Some po- 
sitions (e.g. M31G) suffer from signal subtraction caused by the 
off positions (wobbler mode of observation), and we cannot ex- 
clude that the real signal could be a bit larger. 

5.1 .2. Large-scale tidal streams 

The velocity range measur ed in CO i s comparable with the 
measu rements performed by Ibata et al. (2004); Cha pman et alJ 
(2008) for stellar streams and for extra-pla nar gas and high- 
veloc i ty clouds detected in HI o utside the disc dBraun & Th ilker, 
l2004t IWestmeier etafl I2008I). ICasoli & Combes (1988) have 
detected CO up to 16 kpc. iBraun & Thilkerl (120041) detected a 
faint bridge of HI emission which appears to join the systemic 
velocities of M 3 1 with that of M 33 and continues beyon d M 3 1 
to the north-west. Davies's cloud (We stmeier et all |2008) lies in 
the north-west part of M3 1 's disc and could belong to this bridge. 
Similarly to the stellar streams observed in the south-east part of 



M31's disc dlbata et all 120071) . the HI gas ex hibits large design 
loop- like features (e.g. the Magellanic stream, Kalberla & Haud, 
2006). The second component detected in the molecular gas we 
have detected could be associated to the M3 1 -M33 bridge or any 
gaseous-loop relics. It would then correspond to some accretion 
of gas onto the disc, generating an inner polar ring and lead- 
ing to compression and excitation of the CO. This configuration 
can be compared to th e inner polar gaseous disc discussed by 
ISU'chenkoetalJ d2011l) in NGC 7217, where the disc is face-on. 

While HI studies mention contamination by the Galaxy 
dWestmeier et all [2008), our CO detections are related to ex- 
tinction patterns that are clearly associated to M31. The CO 
detected for M31D has a mean velocity of -74 kms^'is typ- 
ically below the limit for Galactic contamination used for 
HI. However, the typical velocity dispersion of HVC in the 
Galaxy is 8.5 kms" 1 , while the dispersion measured for M31D 
is 24 kms -1 . Furthermore, it exhibits a higher than expected 
Ico/E(B-V) ratio, rather suggesting that it lies behind the bulge 
and not in front of it. Finally, no C O from the Milky Way has 
been detected bv lDame et alJ d200ll) in the vicinity of M3 1 . 

Last, we cannot exclude that some non-circular velocities de- 
tected in CO result from relics of M3 1 's formation: scattered de- 
bris or other extra-planar gas. 

5.2. Scenario 2: a possible bar 

Another a priori plausible scenario could be that the second ve- 
locity component is associated to a possible bar. Bars are fre- 
quent in spiral gal axies (60 to 75% are barred according to the 
bar strength; e.g. IVerlev et all 120071) . and it is likely that M31 
is barred or has bee n barred, because of its triaxial boxy bulge 
(Beat on et all 120071). However, up to now, the bar is not visible 
in the gas in either the morphological or kinematical data. 

5.2.1. Previous bar propositions in M31 

There have been several propositions of bar models for M3 1 . 
However, they are all contradictory (orientation of the bar along 
the minor axis or the major axis), and there is no coherent model 
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that fits the observations. IStark&Binnevld 19941) assume the ex- 
istence of a Ferrers bar potential in the centre to explain some 
non-circular motions there, but there was no complete velocity 
field, and they do not try to fit the observed stellar distribution 
with a bar. Their main point is to propose that the morphology 
of the gas in the centre of M3 1 looks round, in spite of the bar, 
because of the high (77°) inclination of M3 1 : the bar would be 
aligned along the minor axis. The xl orbits, elongated along the 
bar, i.e. the minor axis, are then projected as almost round, be- 
cause their axis ratio could be as high as 4.44= l/cos(77°). The 
authors propose xl orbits between 0.7 and 3 kpc in radius, and 
perpendicular x2 orbits inside 0.7 kpc. The gas would be shocked 
at the transition between xl and x2 orbits, in a ring of 0.7 kpc 
radius. However, as is demonstrated below in Sec. 15.2.21 this 
scenario cannot account for the large observed velocity splitting 
along a single line of sight on the minor axis, at ~ 3.5kpc from 
the centre. 

lAthanassoula & Beaton (2006) investigate the twist of 
isophotes of the NIR 2MASS image, which shows a boxy shape 
bulge. They compare this morphology with simulated stellar 
models with bars, without gas. They conclude that the NIR mor- 
phology of M3 1 in the centre favours a bar, seen almost side-on, 
i.e. the angle between the bar and the major axis of the galaxy 
is between 20 and 30°. The size of the boxy-bulge is 1.4 kpc 
in radius, and the total length of the bar would be 4 kpc in ra- 
dius, because a thin bar is identified beyond the boxy bulge 
dBeaton et all 120071) . 

Let us note that in these two above propositions, the bar al- 
most end-on, or almost side-on, there is no large angle between 
the bar and the symmetry axes of the projected galaxy (major 
or minor axis). This is required because of the morphological 
symmetry of the disc (no bar is obvious in the gas), and also 
because of the fairly regular velocity field. If the bar were in- 
clined by ~ 45° with respect to the major axis, strongly skewed 
(S-shape) isovelocities should be seen in the gas, i.e. the kine- 
matical mino r axis would not be perpen dicular to the major axis, 
as in N2683 dKuzio de Narav et all 120091) . which has been com- 
pared to M31. In the latter, owing to the relative deficiency of 
gas in the centre (iBrinks & Shane! 1 19841 iNieten et all 120061) . it 
is difficult to trace isovelocities, but Figs. [6] and [I0]reveal that 
there is no strong skewness. 

lAthan assoula & Beaton! d2006l) also present a 2D gaseous 
velocity field in an analytical bar potential (not fitted to M31) 
and project it edge-on, i.e. all regions of the galaxy up to 20 kpc 
can thus be seen artificially on the same line of sight, which 
makes a big difference from a 77° projection corresponding to 
M3 1 's real configuration. The model 2D gas morphology clearly 
shows a bar, with emptied regions that are not similar t o the 
de-projected view of M31 gaseous disc (e.g. Braun, 1991). The 
authors compare the position-velocity diagrams parallel to the 
major axis that they derived from the 2D-edge-on model to 
those observed in HI. The latter reveal conspicuous figure-eight 
shapes that turn at the end of the HI extent (R =70'~ 15 kpc), and 
have been interpreted a s caused by the warp dHendersonl 119791 
Brinks & Burton, 1984). These features are much larger than the 
(3-4 kpc) bar extent. Note that in the galaxy NGC 2683, as in- 
clined on the sky plane as M31, the position-velocity diagram 
along the major axis is not a figure-eight shape, but a parallelo- 
gram that ends at a certain radius, precisely the radius of the bar 
R= 2.2kpc= 45" dKuzio de Narav et all 120091) . The situation is 
completely different in M3 1 . 



5.2.2. Special case of M31 

Is it possible to obtain two velocity components on each side of 
the systemic velocity along the same line of sight, near the mi- 
nor axis (position of M31-I and J), at a distance of 210arsec, 
i.e. 3.5 kpc from the centre (if the points are in the main inclined 
plane)? The possibility to gather several components in the same 
beam strongly depends on the spatial resolution. Here, the CO(2- 
1) beam is 12"~ 40 pc on the major axis, however, the inclination 
of the plane limits the resolution. Typically, gaseous discs in the 
centre of galaxies have a characte ristic scale-height of 50 pc, as 
determined in the Milky Way (e.g. lSanders et allll984l) . The in- 
clined line of sight across the disc of M3 1 encompasses a region 
of size =tan(77°) x 50pc= 216pc in the direction parallel to the 
minor axis (while it is the size of the beam, 40 pc, in the perpen- 
dicular direction). The region explored in one line of sight then 
goes along the minor axis from 3.4 kpc to 3.6 kpc, and only 40 pc 
in the perpendicular direction. The only place where it would be 
possible to have two velocity components on each side of the 
systemic velocity is the very centre, because the 216pc region 
explored would encompass the two sides of an elongated orbit. 
Assuming that the bar is inclined with some angle with respect 
to the major axis, the elongated orbits will provide red- and blue- 
shifted velocities on each side of the centre. This is the case in 
NGC2683 for instance, which is revealed by a long slit exactly 
along the major axis. It is possible in the centre, because the 
bar is inclined with respect to the symmetry axes. However, at 
a distance of 3.5 kpc from the centre (±100 pc), it is impossi- 
ble to obtain these two V-components, with the same category 
of orbits, either xl or x2. The only other possibility so far from 
the centre would be to consider a change of orientation in the 
orbits, or a sudden shock, as occurs from xl to x2, or crossing 
an arm. However, the amplitude of the shock, 260 km s -1 , and 
the positions of the two components on the two sides of V is 
unrealistic, especially for a possible weak bar. 

5.2.3. Examples of bar shocks 

Several strongly barred galaxies have been studied in Ha spec- 
troscopy, determining the kinematics with high spatial resolu- 
tion, comparable to what we have in the CO gas towards M3 1 . 
NGC 1365 is the prototype of a very strong bar. The position an- 
gle of the bar does not coincide with the symmetry axis of the 
projection, so that the bar signatures are obvious. With an in- 
clination of i=42°, and a distance of D=20Mpc, the observed 
velocity gradient across the bar dust lanes is determined to be 
at maximum 40kms~ 1 , with a spatial resolution below lOOpc 
dLindblad et all 119961) . If deprojected, the maximum would be 
60 km s~ 1 . NGC 1 300 is another remarkable and strongly barred 
spiral galaxy, where the kinematic s has been studied in detail in 
Ha, on scales smaller than 100 pc. Lindblad et all d 1997b present 
a projected velocity gradient across the dust lanes with shocks 
of 20-30 km s 'projected on the sky. When deprojected, this be- 
comes a maximum of 35-50 kms -1 . To reach these strong ve- 
locity gradients, a highly contrasted gas structure must be found, 
which is not the case for the M3 1 central region. 

In the molecular component, interferometric observations of 
nearby barred galaxies can now yield sub-arcsec resolution, cor- 
responding to scales smaller than lOOpc, th erefore comparable 
to the single-dish resolution towards M3 1 . Garcia-Burillo et all 
d2005l) have presented non-circular motions caused by bars, 
which are of the same orders of magnitude ~ 60-70 km s . 
More recently, w it h eve n higher sensitivity and resolution, 
van der Laan et al .] (1201 lh have for the first time seen dedou- 
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bled velocity profiles in the CO data towards the barred galaxy 
NGC 6951. If interpreted in terms of crowding of orbits, they 
reveal velocity differences of 80 km s _1 (within 500 pc), which 
in deprojection would amount to 110 kms -1 . The CO velocity 
profiles are, however, very broad in the centre of this galaxy, 
showing a starburst ring, and it is not sure what is caused by 
turbulence or by the bar. The molecular gas surface density in 
the centre is on average 100M o /pc 2 (corresponding to a total 
mass of 1.5 x 1O 9 M inside 2kpc), which may explain the un- 
stable and turbulent gas-forming stars, while in M3 1 the surface 
density is at least one order of magnitude lower. The velocity 
dispersion in the distinct velocity components observed in CO is 
quite low (15 kms -1 ), which corresponds to a stable gas layer. 
We conclude that the bar is not a likely explanation for the spe- 
cial kinematical phenomenon observed towards the M3 1 centre. 

5.2.4. Other interpretations 

From the Hi-derived gas morphology, iBraunl (1 1 99 lb suggests 
a two-armed trailing density wave with a pattern speed of 
15 kms 1 kpc~', leading to a corotation at 16kpc, an ILR at 
5 kpc, and an OLR at 22 kpc. He proposes as the driver of the spi- 
ral structure the collision with M32, which could also be respon- 
sible for the many tilts of the M31 plane. The nea rly head-on 
collis ion with M32 has been developed in detail bv iBlock et al.l 
(2006), which allows them to account for the two off-centred gas 
rings observed in M3 1 at 0.7 kpc and 10 kpc radius. Indeed, these 
authors note that the gas morphology in M3 1 does not show the 
usual signatures of the response to a bar. Although it is pos- 
sible that the galaxy disc possesses a stellar bar, aligned with 
the triaxial bulge, oriented roughly along the major axis, the ob- 
served 0.7 kpc dust ring is aligned almost along the minor axis, 
and not aligned with the potential bar along the major axis. The 
size of the inner ring does not correspond to a possible inner 
Lindblad resonance, nor the outer ring at 10 kpc radius to a cor- 
related outer resonance. Moreover, the inner disc appears highly 
inclined with respect to the main disc, and the inner ring is off- 
centred by 40% of its radius, which is not expected in barred 
galaxies. This offset centre and the tilted plane of the inner disc 
strongly point towards a perturbed origin of the inner features of 
M31. 



The resulting configuration is similar to the scenario pro- 
posed in Sect. 15.1.21 Below, we discuss some modelling that 
enables one to explain the observations. We stress that this mod- 
elling is meant to be schematic rather than reproducing the ob- 
servations in a fully self-consistent manner. 

We used a model for the gas, through static simulations of 
dynamical components represented by particles, with different 
geometrical orientations on the sky, but embedded in a gravita- 
tional potential representing the observed rotation curve of M3 1 . 

We represent the gaseous disc of M3 1 by a fairly homoge- 
neous Miyamoto-Nagai disc of particles (with radial scale of 
1 kpc, and height of 0.2 kpc) to be able to vary easily the inclina- 
tion and position angle as a function of radius. We used typically 
half a million particles to have sufficient statistics. We plunged 
the gas disc into a potential made of a stars and a dark matter 
halo. The stellar component is composed of a bulge and a disc. 
The bulge is initially distributed as a Plummer sphere, with a 
potential 

* 4(r) = __^ => (1) 

V r2 + r l 

where Mb and r/, are the mass and characteristic radius of the 
bulge, respectively (see Table|6]i. 

The stellar disc is initially a Kuzmin-Toomre disc of surface 
density 

Z(r) = Z (l+r 2 /^r 3/2 (2) 

with a mass Mj, and characteristic radius rj. 

The dark matter halo is also a Plummer sphere, with mass 
Mqm and characteristic radius ruM- A summary of the adopted 
parameters is given in Table [6] The resulting rotational velocity 
curve is rising until a maximum of 300kms _I , reached at 3 kpc, 
and then slightly decreases, remaining close to 300 kms 1 in 
the central regions of interest here. This correspon ds to the rota- 
tional velocity observed (e.g. Cari gnan et al.Ll2006h . All molecu- 
lar clouds we aim to reproduce are inside the radius of maximum 
rotational velocity. 

The gas particles are distributed in this potential, with the 
velocity dispersion corresponding to a Toomre Q parameter of 
1.3. 



5.3. Scenario 3: head-on collision with M32 

iJacobv et al.l (|l985j) have shown that the inner disc of M31 is 
tilted with respect to the main disc. While the inclination of the 
large-scale disc is 77°, the nuclear disc, of size ~ 1 kpc , is in- 
clined by only about 40° on the sky plane. IBlock et al] ([2006) 
discovered in the Spitzer-IRAC maps an inner dust disc of scale 
lkpc by 1.5 kpc, which also appears with this low inclination. 
The ring morphology of this tilted structure together with the 
ring (lOkpc) morphology of the large-scale disc led these au- 
thors to propose a head-on collision with a small companion, 
with a probable candidate being M32. 

In this scenario, we assume that the peculiar component 
comes from a tilted ring-like material, likely coming from the 
perturbed gas from the recent M32 collision. The collision might 
have perturbed the nuclear disc, which was already tilted, and 
produced some local warp before the annular density wave prop- 
agated outwards. Besides its interest for the inner configuration 
of M3 1 , this scenario also proposes an explanation for the 1 kpc 
ring, resulting from the propagation of the initial annular wave. 



Table 6. Model parameters 



Component 


radius 


Mass 


Mass fraction 




[kpc] 


[M ] 


[%] 


Bulge 


0.2 


1.6el0 


4.5 


Disc 


4.0 


7.el0 


20. 


Halo 


10. 


27.el0 


75.5 



To reproduce the observations, we adopted an inclination of 
43° for the gaseous nuclear disc, inside the 1.5 kpc radius. At 
the boundary of this disc, we assumed a progressive warping of 
the disc plane, so that the inclination on the sky grows from 43 
to 77 degrees, over 300pc. The details of this transition, how- 
ever, are not constrained by the observations, because the lat- 
ter are all inside the 1 kpc radius. As for the position angle on 
the sky, the nuclear disc has PA= 70°, unlike the main disc, 
which has PA= 35°. In projection on the sky, the nuclear gas 
disc model gives velocities that agree with the observations, at 
least with the main velocity component. In many observed points 
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Fig. 12. Mean density of CO emission from the homogeneous 
nuclear disc plus tilted ring of the model. The field of view is 
200 arcsec in radius, or 0.76 kpc in radius. The density scale, 
indicated in the wedge, is in arbitrary units. The positions of our 
CO observations are indicated with circles corresponding to the 
beams. 




100 -100 

ARA (") 



Fig. 14. Density-weighted velocity dispersion (second moment 
of the simulated velocity field). The wedge gives the veloc- 
ity scale in kms . The locations of double components with 
counter-rotating features correspond to the blue and purple re- 
gions. 




Table 7. Geometrical parameters 



Fig. 13. Density-weighted mean velocity (first moment) of the 
simulated nuclear region in the same field of view as Fig fT2l 
Signatures of the tilted ring can be seen at the boarder of the 
field. The wedge gives the velocity scale in kms" 1 . 



only this main component is observed, with a broad line-width 
(FWHM=50-70 kms 1 ), located on the blueshifted side in the 
SW and on the redshifted side in the NE of the major axis. In 
some of the points there is an additional peculiar velocity com- 
ponent, located on the opposite side (redshifted in the SW), and 
narrower (20 kms -1 ). 

Trying several parameters, we found that material on a ring- 
like orbit, with 40° inclination, and PA= -35° gives kinemati- 
cal features compatible with the observations (cf Figure [TSUI 5 1 ). 
The modelled ring has a width of 0.4 kpc, and a mean radius 
of 0.6 kpc, coinciding with the observed dust ring. The adopted 
position angle and inclination of the various components are dis- 
played in Table [7] 



Component 


Main disc 


Nuclear disc 


Inner ring 




[°] 


[°] 


[°] 


PA 


35 


70 


-35 


Inclination 


77 


43 


40 




-100 100 
Velocity (km/s) 



-100 100 
Velocity (km/s) 



Fig. 15. Spectra extracted from the simulated cube at the (RA, 
DEC) offset in arcseconds indicated in each panel. The vertical 
scale is in arbitrary units. The velocities are centred at the sys- 
temic velocity. These offsets correspond the M31G and M31I 
positions. The dumpiness of the gas (not included in the mod- 
elling proposed here) is expected to modulate the relative inten- 
sity of the various components additionally. 



To draw these maps, we built data cubes corresponding to 
the observations, with a pixel of 6 arcsecond, and channels of 
12 kms 1 , and the data smoothed to a beam of 12". The pixel 
size of the cubes are therefore (60, 60, 40) to describe the inner 
parts studied here in CO lines and dust extinction. To take into 
account that the gas and dust distribution is not ho mogeneous, 
but pa tchy, we smoothed the dust map obtained by iBlock et aU 
(2006) from Spitzer-IRAC images (the 8/im map with the 3.6/im 
contribution subtracted) to the same spatial resolution, and used 
it as a multiplicative filter for our homogeneous disc map. This 
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Fig. 16. Schematic view of the interpretation proposed for the 
CO velocities observed. The inner disc is presented with a PA 
of 70deg and an inclination of 43deg. The inner ring is superim- 
posed with a similar inclination but a position angle of -35deg. 
The straight line indicates the position of the major axis of the 
main disc inclined by 77deg with a PA of 35deg. The bars la- 
belled "N" and "E" have a length of 150pc. The near (resp. 
far) sides of the two components are coloured (resp. kept white). 
The red (resp. blue) colours correspond to the redshift (resp. 
blueshift) relative to the systemic velocity. 



takes into account that the dust inner ring is offset by 0.5 kpc 
(131") from the centre of the galaxy. 

Typical spectra were extracted in the zone where double fea- 
tures are expected at the offset indicated in Fig[l5]in arcseconds. 
On the one hand, the broad characteristics of the blue-shifted 
velocity component is caused by beam smearing of the velocity 
gradient, as illustrated here with a 12" beam. On the other hand, 
because the ring is narrow in radius (~ 170pc) and in a region 
of a flat rotation curve, the velocity gradient in the beam is very 
flat and a narrow line is observed in the redshifted component. 

In addition, some emission at the systemic velocity is ex- 
pected because of the inclination of the 77deg main disc. It 
is therefore not necessary to evoke the large-scale warp to ex- 
plain the weak emission at the systemic velocity, which can 
be accounted for by the inner main disc. For positions like 
M31I/M31J/M31G located in the north-west ring, main disc re- 
gions at 3.5 kpc project in the same beam. 

6. Discussion and conclusions 

We have detected CO emission, tracing the presence of molec- 
ular gas in several positions in the central kpc of M31, corre- 
sponding to dust extinction. We have also shown that there is 
an excellent correspondence between the dust extinction mea- 
sured in the optical and the dust emission in the near-infrared 
at 8^m (right panel of Figure ITTb . Some of the observed posi- 
tions correspond to the conspi cuous inner ring, d etected in dust 
emission by IRAC on Spitzer dBlock et"aL . 2006|). This ring was 
also present in ISO data studied by Willaime et al. (200lj) and 



velocity components are detected on each side of the systemic 
velocity. The component with the expected velocity, according 
to the rotation curve and the positive position angle, is broad, and 
is likely to correspond to the nuclear disc. Exploring a wide re- 
gion of the strong central rotational gradient, the beam-averaged 
spectrum may be lOOkm/s broad. The second component, red- 
shifted on the other side of the systemic velocity, has a small 
velocity width, because it is not a full disc, but a relatively nar- 
row ring. We used radii from 0.4 to 0.8 kpc for the model. 

We summarise in Figure [16] the geometry of the two com- 
ponents we propose in our modelling. The peculiar component 
appears to be counter-rotating, because the gas is in a tilted ring, 
almost perpendicular to that of the inner disc. Its inclination on 
the plane of the sky is similar, which explains why the amplitude 
of the projected rotation is the same as for the regular compo- 
nent. In the figure the arrows indicate the apparent sense of rota- 
tion. The two discs appear to rotate in the same direction, which 
is also that of the main M3 1 disc. The apparent counter-rotation 
appears only in the NW (and SE) regions, where blue and red re- 
gions superpose. The winding sense of the spiral structure in the 
observations is also sketched. It is possible to see that the arms 
are trailing, exactly like the M31 main disc. 

In summary, the apparent counter-rotation is only caused by 
a warping and distortion of the central components, possibly 
triggered the proposed head-on collision with M32, but is not 
a true counter-rotation. 
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Fig.A.l. Signal detected in 2000 at the J2000 coordinates RA 
00:43: 18.5 DEC +41 : 17:06, in the OFF signal of position switch 
observations of M31C. This spectrum, corresponding to an ex- 
posure time of 80 min, was reduced to a 24" beam. The y axis 
displays antenna temperature T^. 
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Fig. A.2. Signal detected in the OFF positions of M31D obser- 
vations, with the wobbler switch observing mode, in 1999 (wob- 
bler throw: 157"- 165"). The y axis displays antenna tempera- 
ture T* . 



Appendix A: CO signal detected in the OFF 
positions 

We provide here a short description of the CO signal detected 
in the OFF positions of the IRAM observations presented in this 
paper. Their interpretation is beyond the scope of this paper, be- 
cause they explore a larger area than the ON positions and are 
only partial. 



OFF position of M31C (position switch) 

A clear signal was detected at Vb = -402.3 km s _1 (cr = 
12.5kms _1 , I co = 0.85Kkms"') and V = -320kms"' (cr = 
5.7kms~', Ico = 0.36Kkms~') in the OFF position of the po- 
sition switch observations of M31C, as displayed in Figure |A~T1 
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Fig. A.3. Signal detected in the OFF positions of M31G obser- 
vations, with the wobbler switch observing mode, in 1999 (top: 
wobbler throw 75"- 90") and in 2000 (bottom: wobbler throw 
94" - 104"). The 2000 observations (bottom panel) was reduced 
to a 24" beam. The y axis displays antenna temperature T* A . 



OFF position of M31D (wobbler switch) 

For M31D no OFF signal was detected in the position switch 
scans performed at the OFF position 122", 321". However, there 
is a detectable OFF signal for the scans obtained in wobbler 
switch mode in 1999, as displayed in Figure IA~2l This signal 
corresponds to a wobbler throw centred with 157" and 165" 
from the ON position. Because the observing procedure was de- 
signed to avoid extincted areas and the area surrounding M3 1 D 
was not crowded, we suspect that this OFF signal close in ve- 
locity to the ON signal is caused by gas on the back side (not 
detected in extinction). However, complementary observations 
are required to confirm this hypothesis. The detected signal is at 
y = -^^kms" 1 (cr = 31.4210ns- 1 , l co = l.^Kkms" 1 ) 
with an exposure time of 162 min. 

OFF position of M31G (wobbler switch) 

No OFF signal was detected in the position switch scans ob- 
served for M3 1G in 2000 and only wobbler mode has been used 
in 1999. However, a signal was observed in the OFF positions 
obtained with the wobbler switch mode, as displayed in Figure 
IA.3I The throw positions were different in 1999 and 2000 and 
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the OFF signal detected was also different. In 1999 the detected 
signals are at Vo = -204.5 km s _1 (cr = 17.5kms _1 , Iqo 
2.07 K km s ) and Vo = -242.610ns- 1 (cr = 23.2 km s" 1 , I C0 = 
1.42Kkms _1 ) corresponding to an exposure time of 180min. 
In 2000 the OFF signal is detected at V = -303.3kms _1 
(cr = 17.6kms _I , Ico — 1.36Kkms~ 1 ), corresponding to an 
exposure time of 142 min. As discussed in the paper, these OFF 
positions superimpose on a wide ON signal. Because we do not 
know the position of the wobble throw (only its distance from 
the ON position), we can only argue that these OFF signals de- 
tected within 1.8 arcmin from the ON position can be associated 
to the same structure. 

Appendix B: HI spectra 

Figure |B~T1 displays the HI velocities obtained for the different 
positions observed in CO. The positions M31I and M31J exhibit 
one possible component compatible with the disc component de- 
tected in CO. There is no clear HI counter-part for the ring com- 
ponent. 
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Fig. B.l. HI spectra obtained with a 1 arcmin resolution for the CO positions discussed in this paper. In each spectrum the velocity 
of the CO emission is indicated in red. 
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